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I investigate the decoherence of two-mode squeezed vacuum states by analyzing the relative 
entropy of entanglement. I consider two sources of decoherence: (i) the phase damping and (ii) 
the amplitude damping due to the coupling to the thermal environment. In particular, I give the 
exact value of the relative entropy of entanglement for the phase damping model. For the amplitude 
damping model, I give an upper bound for the relative entropy of entanglement, which turns out to 
be a good approximation for the entanglement measure in usual experimental situations. 
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Quantum entanglement is an essential ingredient in quantum communication and computation jjj. Therefore, it is 
of great importance to quantify the entanglement to assess the efficacy of quantum information processing. Recently, 
unconditional quantum teleportation of an unknown coherent state has been realized experimentally by exploiting a 
two- mode squeezed vacuum state as an entanglement resource Q shortly after the theoretical proposal of Ref. ||. 
_ -, i The two-mode squeezed vacuum state shared between two parties - Alice and Bob - is formally generated from the 
vacuum state \vac) by the unitary transformation, U{r) = exp [-r (a\4 - &1 a 2 )] , where r(> 0) is called a squeezing 

parameter. The indices 1 and 2 refer to the optical mode of Alice and that of Bob, respectively. The two-mode 
squeezed vacuum state \^f) = U(r) \vac) is written in the Fock basis as \^) = (coshr) 1 ^^L tanh™ r |n, n), where 
\n, n) = \n) 1 ® \n) 2 . Since is a pure state, its amount of entanglement is uniquely quantified by the von Neumann 
• entropy of the reduced state of Alice or Bob: 
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8 : 5(|^)) = cosh^log 2 (cosh^)-sinhVlo g2 (sinhV). (1) 



In real experimental situations, due to coupling to the environment, the entangled state inevitably loses its purity; it 
becomes mixed. This phenomenon - decoherence - is the most dangerous obstacle for all entanglement manipulations. 
G \ Several protocols for entanglement enhancement or purification in continuous variable systems have been proposed 
^p H ' H~H an d the decoherence of continuous variable states has been studied by making use of Bell's inequality [Q. 

However, the quantification of entanglement of mixed entangled states in continuous variable systems is still not well 
^ ■ understood. 

In this paper, I investigate the decoherence of two-mode squeezed vacuum states by analyzing the relative entropy 
of entanglement . I consider two sources of decoherence separately; (i) the phase damping and (ii) the amplitude 
| damping due to the coupling to the thermal environment. The damping is assumed to affect each mode of the state 
independently with the same coupling parameters. In particular, I show that the exact calculation of the relative 
entropy of entanglement can be performed for the phase damping model. For the amplitude damping model, I 
give an upper bound for the relative entropy of entanglement, which turns out to be a good approximation for the 
entanglement measure in usual experimental situations. 

The relative entropy of entanglement of pure states reduces to the von Neumann entropy of the reduced state of 
either subsystem. For a mixed state p it is defined as Er(p) = muio-gp S^pHcr), where 5 f (p||cr) = Tr [p (log 2 p — log 2 a)] 
is the quantum relative entropy. The minimum is taken over T>, the set of all disentangled states. It is usually difficult 
to calculate the relative entropy of entanglement for mixed states, except for some specific states. Recently, the 
following theorem on the relative entropy of entanglement has been proved |9) and it turns out to be quite suitable 
for my analysis. 

Theorem 1 For a bipartite quantum state described by the density matrix of the form, 

P = ^2 a »i.«2 l^nu^m) (0n 2 !^n 2 | , (2) 

111 ,™2 

the relative entropy of entanglement is given by 

Er(p) = a„,„ log 2 a n , n + Tr (p log 2 p) , (3) 
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and the disentangled state p* which minimizes the quantum relative entropy S (p\\p*) is 



P* = y^an.n \4>n,1pn) (<f>n,1pn\ 



(4) 



Here, \<f) n ) and \ij) n ) are orthonormal states of each subsystem. 

First, I consider the phase damping model. The density matrix obeys the following master equation in the interaction 
picture, 



dt 



P (t) = (C 1 + C 2 )p(t), 



with 



r T 
C lP =- 



2a\a,ipa\ai - (a\a^j p- p (a\a^j 



The solution of Eq. (||) with the initial condition p(t = 0) = l^) (\&| is calculated as 



Pit) 



cosh 



oo 

(tanhr)™ 1+ ™ 2 exp (—jt \ni - n 2 \ 2 ^j |ni,ni) (n 2 , n 2 \ 



(5) 



(6) 



(7) 



It should be noted that the density matrix of Eq. (Q) takes the form of Eq. (|^) in Theorem [l]. Consequently, it is 
possible to calculate exactly the relative entropy of entanglement of the state of Eq. (Q) . Figure |l| shows the relative 
entropy of entanglement Er thus computed as a function of squeezing parameter r and the degree of damping d = jt. 
In numerical computations, the truncated photon number has been taken to be max(ni) = max (n 2 ) = 100, the 
value of which is sufficiently large for numerical convergence. It is seen that with an increasing amount of squeezing 
the amount of entanglement decreases rapidly with the damping. For large values of d, Er is vanishingly small but 
remains finite; it is still larger than the conceivable numerical errors. It is not clear from the present numerical analysis 
whether the state is always entangled for finite "ft. 




FIG. 1. The relative entropy of entanglement of the state of Eq. (J7j) as a function of squeezing parameter r and the degree 
of damping d(— "ft). 

Next, I consider the amplitude damping model. The master equation for the density matrix is the same as Eq. (@), 
but 

A'P = Z (1 + n) \ 2a,ipa\ - a\a t p - pa] a^j + \ 2a\pai ~ aia\p - P^afj , (8) 

where n is the average photon number of the thermal environment. To solve the master equation, I first represent the 
two- mode squeezed vacuum state as a continuous superposition of two- mode coherent states M : 
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l*>-/*«"»0.'>fc«0. (9) 

where G(a,r) = exp [— (cothr — 1) |a| 2 ] / (7rsinhr). According to the treatment of Ref. |lC| ], the solution for the 
density matrix with the initial condition p(t = 0) = |^) (^j is calculated as 



p(t) = J da J d 2 pG(a, r)G{J3, r)ai (a, (3) a 2 (a*,/3*) , (10) 

where 

0) = exp Dt (6) exp («ajai) D\ (&) . (11) 



In Eq. @, = n(l - e^ 7 *), A (0) = ex P Ki a ! ~ €j a i) > and « = In {n(t) / [n(t) + 1]} . Furthermore, fi and 

f 2 are given by £1 = e"^ {[n(t) + 1] a + n(t)/3} / [2n(t) + 1], and = e - ^ {n(t)a + + 1] /3} / [2n(t) + 1]. By 
noting the completeness of coherent states, it is straightforward to show that 



i (mi| A (&) exp (Kafai) dJ (6) |m 2 ) i = -^J r 

V / 7T Z V Til'? 

x exp 

Using Eq. @, we obtain 



5 (l7i| 2 + |72| 2 + |7i " Cil 2 + 172 ~ 6I 2 ) + (71 - £1)* (72 - 6) e K + V^Tlm (frtf) - V^Tlm (frtf) 



(12) 



(ni,«3|p(t)|n2,Tu) = 4 ^ 4 \ , } , , / / d 2 72 / As f d 2 ^ (j^ if {llT* 
7r 4 smh r V m!n 2 !n3!n4! J J J J 

x exp [-I71I 2 - | 72 | 2 - | 73 | 2 - | 74 | 2 ] exp [P (7*7* + 7274) + Q (71*72 + 7 3 *74)] , (13) 

where 

P = Re- 1 * cothr, (14) 

Q = {n(t) + Re-** [n(t) + 1 - e"^] } , (15) 

n(t) + 1 

and 

R = |coth 2 r[n(t) + l] 2 - [n(t) + 1 - e~ 7 *] 2 } . (16) 

The integral of Eq. (|l^) can be evaluated by expanding exp [P (7173 + 7274) + Q (7*72 + 7374)] with respect to 7i*7|, 
7274, 7*72, and 73*74. Finally we obtain 

00 

P(t)= C nln 2 K> n l> ("2,n 2 | 

Til ,^2=0 

00 oc 00 00 

+ J2 c n 1 ) ,r l2 l«i,«i+fc)(«2,n 2 + fc|+^ J2 C S«2 \ni + k,m) (n 2 + k,n 2 \, (17) 

k—1 Til ,712—0 fc= 1 711 ,^2—0 

where 

Jfc) _ V"i'"2!(ni + fc)!(n 2 + fc)! ni+n k mm J^'" 2) 1 / Q\ 2 ' 

ni '" 2 sinh 2 r & I!(J + *)!(ni-i)!(na-0! VP/ ' 

To estimate the relative entropy of entanglement of the state pit) of Eq. ([l7|), I write /?(£) as a convex combination of 
new density matrices, po(t) and p^\t) (k = 1,2, ...): 
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pit) = Po po(t) + J^p k +) p^ (t) + f>jf Vjf >(*)• (19) 



k=l k=l 



The first, second, and third terms of the right-hand side of Eq. ( |l9| ) correspond to the first, second, and third terms 
of the right-hand side of Eq. (p7|), respectively. By convexity of the relative entropy of entanglement ||, we have that 

oo oo 

Er (p(t)) < p E R (p Q {t)) + Y,pI +)e * (pI^W) + X>jf ^* • ( 2 °) 

k=l k=l 

The density matrices po(t) and (t) (k = 1,2, ...) have the form of Eq. (||) in Theorem |l|, so the right-hand side 
of Eq. @) can be calculated exactly and it yields an upper bound for the relative entropy of entanglement of p(t). 
At jt = 0, this upper bound E* R gives the exact value of the relative entropy of entanglement of Eq. (|l|). Figures 
H and Hshow E R thus computed as a function of squeezing parameter r and the degree of damping d(= jt). In 
Fig. H (||), n, the averaged number of photons in the thermal environment, has been taken to be 0.01 (0.1). The 
dotted line on the data surface in each figure indicates the separability-inseparability border line, which is given by 
the necessary and sufficient separability criterion for the two-mode squeezed state in the thermal environment jTl| , [l2| ; 
7< = In [l + (l — e~ 2r ) / (2n)] . Within the inside region (above the dotted line) the state is entangled, while within 
the outside region (below the dotted line) it is separable and the relative entropy of entanglement Er vanishes. It is 
seen that the values of the upper bound E* R on these border lines for r < 1 are already negligibly small. Although 
E R is merely an upper bound for Er, it may be considered as an approximation for Er. Since r < 1 and n <ti 1 in 
usual experimental situations, this approximation is a fairly good one in the sense that the upper bound E R « for 
the separable region. 

In summary, the decoherence of two-mode squeezed vacuum states has been investigated by analyzing the relative 
entropy of entanglement for (i) the phase damping model and (ii) the amplitude damping model. For the phase 
damping model, a method for the exact numerical computation of the relative entropy of entanglement has been 
established. For the amplitude damping model, a good approximation for the relative entropy of entanglement in 
usual experimental situations has been introduced. 




FIG. 2. The upper bound for the relative entropy of entanglement of the state of Eq. rt!7) as a function of squeezing 



parameter r and the degree of damping d(= -yt). The averaged number of photons in the therma 
dotted line is the separability-inseparability border line. 



environment n = 0.01. The 
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FIG. 3. The same as Fig. |, but n = 0.1. 
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